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O I. INTRODUCTION 

O ; 

Exact solutions in quantum mechanics are hard to come by. This is especially so in the case of relativistic wave 
equations. In particular, for the Dirac equation only a few exactly solvable electromagnetic field configurations such 
as homogeneous magnetic fields lj, homogeneous electrostatic fields Q, constant parallel magnetic fields etc. are 
known. 

While exact solvability is desirable, in practice it is not always possible to determine the whole spectrum. Recently, 
■ in non relativistic quantum mechanics a new class of potentials which are intermediate to exactly solvable ones and 
non solvable ones have been found . These are called quasi-exactly solvable (QES] problems for which it is possible 
t-H , to determine algebraically a part of the spectrum but not the whole spectrum 0, 0, 0, @ • Usually a QES problem 
admits a certain underlying Lie algebraic symmetry which is responsible for the quasi-exact solutions. Previous 
studies in QES problems deal mainly with the non-relativistic equations, and only rather recently did discussions of 
the quasi-exactly solvability of the Dirac equations appear. Particularly, it had been shown that the Dirac equation 
with Coulomb interaction supplemented by a linear radial Lorentz scalar potential 0, and planar Dirac equation 
with Coulomb and homogeneous magnetic fields |ld| | are QES systems. The Pauli equation and the Dirac equation 
I/"") ' coupled minimally to a stationary vector potential were also shown to be QES More recently, added to the list 
are the Dirac-Pauli equation coupled non-minimally to external electric fields [lj , and Dirac oscillator with Coulomb 
' interaction supplemented by a linear radial Lorentz scalar potential 13]. 

Interest in the Dirac equation with Lorentz scalar potential was mainly motivated by the MIT bag model of quark 
confinement 0, 0] . Unfortunately, QES and exactly solvable Lorentz scalar potentials in such Dirac equation are 
rather scanty in the literature. The system considered in provides the first example of a QES system. Motivated 
_C • by this work, in pl| a reparametrization of the radial Dirac equation with a scalar and a Lorentz scalar potential was 
used to show the existence of infinitely many potentials admitting two QES states. However, explicit construction 
of the QES potentials was not given, as the reparametrized equations are difficult to handle in practice. Later, a 
subclass of QES potentials admitting doublets and multiplets were constructed by restricting the forms of the wave 
functions 1171 . 

In [111 Il2| we have developed a procedure to construct exactly-solvable potentials and s^(2)-based QES potentials 
admitting any finite number of QES states. The method relies mainly on the factorizability, or equivalently, on 
the supersymmctric structure of the Hamiltonian of the system. In the present paper we would like to extend the 
procedure to the Dirac equation with Lorentz scalar potential. We will determine possible forms of the potential 
allowing factorization of the Dirac equation. Once factorization is achieved, the forms of the exactly solvable and 
QES potentials can be easily determined. 

The organization of the paper is as follows. In Sect. II we outline the procedure presented in [Til [T3| . Sect. Ill 
discusses the exact and quasi-exact solvability of the Dirac equation in one spatial dimension. Sect. IV is devoted to 
the three and four dimensional Dirac equations in Cartesian coordinates, with the potential depending only one spatial 
dimension. In Sect. V the Dirac equation with a Lorentz scalar and a magnetic monopole potential is considered. 
Sect. VI concludes the paper. 
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II. OUTLINE OF THE PROCEDURE 

In this section we shall outline the general procedure we adopted in jllL Il2j | for determining exact /quasi-exact 
solvability of the Dirac and the Pauli equation. This procedure makes full use of the close connection between 
quasi-exactly solvable systems and supersymmetry, or factorizability of the system. 
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Suppose for a Dirac system one can reduce the corresponding multi-component equations to a set of one-variable 
equations possessing one-dimensional supersymmetry after separating the variables in a suitable coordinate system. 
Typically the set of equations takes the form 

|- + W(r)j V^- = £+'0+ , 

~+W{r)\il>+ = £-i/>- , (1) 

where r is the basic variable, e.g. the radial coordinate, and 4>± are ) sa y) the two components of the radial part of the 
Dirac wave function. The superpotential W is related to the external field configuration. may involve the energy, 
the mass of the particle, and perhaps some conserved quantum numbers. We can rewrite this set of equations as 

A~A + ip_ = ei>_ , (2) 
A+A-1P+ = , (3) 



with 



Explicitly, the above equations read 



A ± = ±4- + W, e = £+£- . (4) 
ar 



-Jj+W a TW / )^=c^F- (5) 

Here and below the prime indicates differentiation with respect to the basic variable. Eq. iJSJ clearly exhibits the 
supersymmetric structure of the system. The operators acting on ip± in Eq. JSJ) are said to be factorizable, i.e. as 
products of A~ and A + . The ground state, with e = 0, is given by one of the following two sets of equations: 

A+^\r) = , </4 0) (r) = 0; (6) 
A-ipfir) = , ^i 0) (r) = , (7) 

depending on which solution is normalizable. 

One can determine the forms of the external field that admit exact solutions of the problem by comparing the forms 
of the superpotential W with those listed in Table 4.1 of [Isj . 

Similarly, from Turbiner's classification of the sl(2) QES systems [Bj], one can determine the forms of W, and hence 
the forms of external fields admitting QES solutions based on sl(2) algebra. The main ideas of the procedure are 
outlined below. We shall concentrate only on solution of the upper component ?/>_, which is assumed to have a 
normalizable zero energy state. 

Eq. (JSJ shows that satisfies the Schrodinger equation H-ip- = eip-, with 



H- = A-A^ 
_<P_ 
dr 2 

with 



V(r) , (8) 



V(r) = W(r) 2 - W'(r) . (9) 

We shall look for V(r) such that the system is QES. According to the theory of QES models, one first makes an 
"imaginary gauge transformation" on the function ?/>_ 

V>_(r) = 4>{r)e~ 9{r) , (10) 

where g(r) is called the gauge function. The function <f>(r) satisfies 

d 2 (j>{r) ,d(j)(r) f „ , 2 -, 
+ 2 9 + r( r ) + 9 ~ 9 \ 4>K r ) = <WJ • (11) 
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For physical systems which we are interested in, the phase factor exp(— g(r)) is responsible for the asymptotic behaviors 
of the wave function so as to ensure normalizability. The function <p(r) satisfies a Schrodinger equation with a gauge 
transformed Hamiltonian 

H G = --^ + 2W (r)-^ + [V(r) + W^-W 2 ] , (12) 

where Wo(r) — g'(r). 

Now if V(r) is such that the quantal system is QES, that means the gauge transformed Hamiltonian Hq can be 
written as a quadratic combination of the generators J a of some Lie algebra with a finite dimensional representation. 
Within this finite dimensional Hilbcrt space the Hamiltonian Hq can be diagonalized, and therefore a finite number 
of eigenstates are solvable. For one-dimensional QES systems the most general Lie algebra is si (2) . Hence if Eq. I|12|) 
is QES then it can be expressed as 

H G = J2 CabJ a J b + Caja + constant > ( 13 ) 
where C a b, C a are constant coefficients, and the J a are the generators of the Lie algebra sl(2) given by 



d 

dz 



J + = z 2 —-Nz 



J- = ±. (14) 
dz 

Here the variables r and z are related by z — h(r), where h(-) is some (explicit or implicit) function . The value 
j = N/2 is called the weight of the differential representation of sl(2) algebra, and N is the degree of the eigenfunctions 
(f>, which are polynomials in a (N + l)-dimensional Hilbert space with the basis (1, z, z 2 , . . . , z N ) : 

cf)= (z - z x ){z - z 2 ) ■ ■ ■ {z - z N ) . (15) 

The requirement in Eq. I|13|) fixes V(r) and Wo(r), and Hq will have an algebraic sector with N + 1 eigenvalues and 
eigenfunctions. For definiteness, we shall denote the potential V admitting N + 1 QES states by Vm- From Eqs. i|l(J|) 
and l|15fl , any one of the N + 1 functions ip- in this sector has the general form 



ip- = {z - zi)(z - z 2 ) ■ ■ ■ (z - ZAr)exp I- J W {r)drJ , (16) 

where Zi (i = 1,2, ... , N) are N parameters that can be determined by plugging Eq. I|16fl into Eq. Ijlip, The algebraic 
equations so obtained are called the Bethe ansatz equations corresponding to the QES problem. 
Now comes a crucial observation in our procedure: one can rewrite Eq. i|ltj|) as 



^-=cxp(^-y W N (r,{ Zi })dr\ , (17) 

with 

W N (r, {z^ = W (r) ^ — ^L- . (18) 
i=i z% 

There are N + 1 possible functions Wwir, {zi}) for the N + 1 sets of eigenfunctions <f>. Inserting Eq. (|17|l into 
H_ip_ — eip-, one sees that Wn satisfies the Ricatti equation 

Wlf - W' N = V N - e N , (19) 

where is the energy parameter corresponding to the eigenfunction ^_ given in Eq. I|l(j|) for a particular set of N 
parameters {z{\. 

From Eqs. (JSJ, JHJ and 119() it is clear how one should proceed to determine the external fields so that the Dirac 
equation becomes QES based on sl(2): one needs only to determine the superpotentials W(r) according to Eq. i|19|) 
from the QES potentials V(r) classified in p|. This is easily done by observing that the superpotential Wo corre- 
sponding to N — is related to the gauge function g(r) associated with a particular class of QES potential V(r) 
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by 9'( r ) — Wq(t). This superpotential gives the field configuration that allows the weight zero (j — N = 0) state, 
i.e. the ground state, to be known in that class. The more interesting task is to obtain higher weight states (i.e. 
j > 0), which will include excited states. For weight j (N = 2j) states, this is achieved by forming the superpotential 
Wff(r, {zi}) according to Eq. i|18|) . Of the N + 1 possible sets of solutions of the Bethe ansatz equations, the set of 
roots {zi, Z2, ■ ■ . , Zfj} to be used in Eq. (|18|l is chosen to be the set for which the energy parameter of the corresponding 
state is the lowest. 

In the following three sections, we will employ this procedure to classify and construct QES Lorentz scalar potnetials. 



III. TWO-DIMENSIONAL DIRAC EQUATION IN CARTESIAN COORDINATES 

Let us consider a (1 + l)-dimensional Dirac Hamiltonian of the kind 

H = ap + f3(m + V s (x)) , (20) 

where m is the mass of the fermion, p = —id/dx, a and (3 are the Dirac matrices, and V s is the Lorentz scalar 
potential. Such model is of interest in the theory of nuclear shell model [THl ITflj , and as a model of the self-compatible 
field of a quark system |20j . Furthermore, if there exists zero energy solution, then the theory possesses spectral 
asymmetry and fractional fermion number |15ll2l| . 

It has long been known that the system is supersymmetric [22l l2^|. This can be easily shown as follows. Let the 
Dirac matrices be represented in terms of the Pauli matrices as 

a x = o- 2 , P = £7i . (21) 
In this representation the Dirac equation Hip = Etp for the 2-component wave function 

«»>-(&*>) (22) 

takes the form: 

Tx + u{x) ) *~ = m+ ' 

~+U{x)^+ - . (23) 

Here U(x) = V s (x) + m and e = E 2 . Eq. I|23|l is now in the factorized form, with U(x) playing the role of the 
superpotential. As such, this system can be dealt with according to the general procedure outlined in the last section. 



A. Exactly solvable cases 

The exactly solvable cases have been classified in Table 4.1 of [1S|. For the present system, there are six types 
of exactly solvable field configurations, namely, (i) shifted oscillator-like; (ii) Morse-potential; (iii) Rosen-Morse II 
(hyperbolic); (iv) Scarf II (hyperbolic); (v) Rosen-Morsel (trigonometric); and (vi) Scarf I (trigonometric) . 

We note here that the case of linear (shifted oscillator) potential has been discussed in [24J , and the Morse case in 





B. Quasi-exactly solvable cases 



In this case, seven classes of QES potential U (x) can be constructed. These classes correspond to Class I to Class VI, 
and Class X in Turbiner's classification 0. We shall illustrate the construction of Class I potentials below. Potentials 
in the other classes can be constructed accordingly. 

The QES potential belonging to Class I has the form 

V N {x) = a 2 e~ 2ax - a [a(2N + 1) + 26] e' 01 * + c(2b- a) e ax + c 2 e 2ax + b 2 - 2ac . (24) 
Without loss of generality, we assume a, a, c > for definiteness. The corresponding gauge function g(x) is given by 

g (x) = -e~ ax + -e ax + bx . (25) 
a a 
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One should always keep in mind that the parameters selected must ensure convergence of the function exp(— g(x)) in 
order to guarantee normalizability of the wave function . The potential V(x) that gives the ground state, with energy 
parameter e = E 2 = 0, is generated by 



with 



V = U 2 - U' 



U (x) = g'{x) = -ae- ax + ce ax + b 



(26) 



(27) 



To obtain the potentials Un(x) which admit solvable states with higher weights j, we must first derive the Bcthe 
ansatz equations. To this end, let us perform the change of variable z — h(x) = exp(— ax). Eq. then becomes 



+ [2az 2 - (2b + a)z -2c] — + 



-2aNz 



(z)=0 



(28) 



For N > 0, there are N + 1 solutions which include excited states. Assuming <fi(z) — Y[iLi( z ~ z i) m Eq. (|28|l . one 
obtains the Bethe ansatz equations which determine the roots Zi's 



z 

2azf - (2b + a)z l - 2c - 2a > — = 0, i = 1, . . . , N , 



Zi ~ Z[ 



(29) 



and the equation which gives the energy parameter in terms of the roots Zi 's 

N 

Zi 



N 1 

= 2 " c E-- 



(30) 



i=l 



Each set of {z{\ determines a QES energy E with the corresponding polynomial <f>. 

As an example, let us construct U\(x) which admits two solutions. This corresponds to the case with N = 1 and 
<f>(z) = z — Z\. According to Eq. (|29|l . the root z\ satisfies 



2az\ - (2b + a)zi - 2c = , 



which gives two solutions 



(26 + a) ± v / (2b + a) 2 + 16ac 



4a 



(31) 



(32) 



The corresponding energy parameters are 



2ac- 



z\ 



(2b + a) T V(2b + a) 2 + 16ac 



(33) 



For the parameters assumed here, the solution with root z 1 = — \z 1 \ < gives the ground state, while that with root 
> gives the first excited state. The superpotential is constructed according to Eq. (|18fl as 



Ux(x) = U 



h'(x) 



h(x) — z 1 



b . 



Zi e" 



(34) 



The ground state and the excited state have energy E 2 = tT — e~ = and E 2 = e + — e~ = a^J (2b - 
respectively. 

Potential Un (x) admitting N + 1 solvable states can be constructed accordingly. 



16ac, 
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IV. THREE- AND FOUR-DIMENSIONAL DIRAC EQUATIONS IN CARTESIAN COORDINATES 

A. Three-dimensional Dirac equation 

A (2 + l)-dimensional Dirac Hamiltonian has the form 

H = a ■ p +/3(m + V s (x)) , (35) 

where p = (j> x ,p v ) = —i(d/dx,d/dy), and a, (3 are Dirac matrices. The potential V s (x) is assumed to depend only 
on x. 

The Dirac equation can also be cast into factorized, or supersymmetric form. We represent the Dirac matrices in 
terms of the Pauli matrices <r as follows: 

a x = <t x , a y = <Jy , (3 = a z . (36) 
As V(x) depends on x only, the wave function can be written as 

^ = elkyV (uU) ' (37) 

where k y is a real constant, and f± are real functions of x. The Dirac equation becomes 

U(x) p x ~ik y \ ( f_{x) \ _ _ / f-(x) \ , . 

Px + t k y -u{ X ) ) \ u(x) j-^y u(x) ) ■ ^> 

Here U{x) is again defined as U(x) = V s (x) + m. To cast Eq. I|38|l into supersymmetric form, we transform the wave 
function by a unitary matrix T: 

Ux) \ ( i$_(x) \ = t / f-(x) \ J_ ( 1 i 



f + (x) J '\^ + {x))- T {f + (x))> T -V2^1 ) - (39) 



Then the Hamiltonian transforms as 



"^'"^(-.(-iu^'^x^)- « 4o) 

and the Dirac equation becomes 

■^+U(x)jil>- = (E + k y )^+ , (41) 
— + f/(j;)j^ = (£-y^, (42) 

which is exactly in the same form as Eq. (JTJ . Hence the exact / quasi-exact solvability of Eq. (|42|l can be discussed as 
in the one-dimensional case. 

B. Four-dimensional Dirac equation 

Let us now consider exact/quasi-exact solvability of four-dimensional Dirac equation with a Lorentz scalar potential 
dependent only on one spatial variable, say z. Such system is useful in describing ^-dependent valence- and conduction- 
band edge of semiconductors near certain points in the Brillouin zone |2q . 

The Hamiltonian has the form 



H = a-p + P(m + V s (z)) 

a p — ill(x) 

a ■ p + iU(x) 



(43) 



where p = (p x ,p y ,p z ) — —i(d/dx, d/dy, d/dz ). Here U{z) = V s (z) + m as before, and we have chosen the Dirac 
matrices in the supersymmetric representation 26] 
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To see that the Hamiltonian is factorizable, it is more convenient to consider its square H 2 : 

H~ 
H 

where 

HT = -W 2 + U 2 (z)±U'(z)a z . (46) 

Let the wave function be in the form 

,/, _ ik x x+ik y y ( lp-( z )X- \ ( An\ 

Here k x and k y are real constants, i^±(z) are real functions of z, and \± are two-component eigen-spinors of a z : 
a z x± = ±X±- Then the eigenvalues problem of H 2 ip — E 2 ip reads 

d 2 



, :2 +U 2 TU'j ^=e^, (48) 

with e — E 2 — k 2 — k 2 . Now it is obvious that Eq. (|48() is in the supersymmetric form Eq. (JTJ, and construction of 
QES potentials U(z) proceeds as in the previous cases. 

V. DIRAC EQUATION WITH A LORENTZ SCALAR AND A MAGNETIC MONOPOLE POTENTIAL 

Dirac equation with a spherical Lorentz scalar potential was originally motivated by the MIT bag model intended 
for describing quark confinement [bH Il5| . Now we would like to determine the forms of spherical Lorentz scalar 
potentials which admit exact /quasi-exact solutions. 

Generally, in the presence of an electromagnetic field (V, A) and a Lorentz scalar potential V s , the Dirac equation 
is given by 

H = a ■ (p - ieA) +/3U + eV , (49) 

where e is the charge of the fermion and U = m + V s . Quasi-exact solvability of the system when V = V s — has been 
demonstrated in connection with the Pauli equation [ll|. In the absence of the vector potential A and V, the system 
has been shown in the last two sections to be QES in one spatial dimension. In other coordinates, it is rather difficult, 
if not impossible, to find QES potentials. In this section, we shall consider the system in spherical coordinates. 
When A = 0, the Dirac equation in spherical coordinates can be separated into a two-component radial equation 

153 



l- r + T eV-E-U\ ( f(r) 
eV~E + U + * J \g(r) 



(50) 



The k in the centrifugal term is related to the total angular momentum quantum number j — 1/2, 3/2,... as 

k = ±(J + 1/2). It follows that k is nonzero integer: n = ±1, ± 2, As mentioned in the Introduction, QES 

potentials with both V and U arc difficult to find. The first example was presented in 0, in which V is taken to 
be Coulomb- like and V s is linear in r. Existence of infinitely many QES V and V s with two solutions based on some 
reparametrization scheme was shown in [l(j, and explicit construction of a subclass of such doublet (and multiplct) 
QES potentials were presented in 

We want to determine possible forms of exactly solvable and QES potentials based on factorizability of the system. 
Unfortunately, the Dirac equation Eq. I|5U|I is difficult to factorize, owing to the presence of the non-vanishing cen- 
trifugal term re/r p3l |. Nevertheless, it has been shown in |28| that k could be zero in the presence of a Dirac magnetic 
monopole described by 

I — cos 9 , . 

Ar = A e = , A 4 ,= g — - . 51) 

r sin 

Here g is the magnetic charge of the Dirac monopole. In it is shown that 

— / 2' I' • • ■ if <? = ; 

3 \ \q\-h, kl + l, kl -§,••• if q^O, 1 ' 
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where q = eg — 0, ± 1/2, 1, . . . according to Dirac's quantization condition, and 



X X2 



«=±y(i+§j -? 2 - ( 53 ) 

Hence, k = when j = |g| — 1/2, which is possible only when the magnetic charge g / 0. In this case, the Dirac 
equation with a scalar and a Lorentz scalar potential in the presence of the Dirac monopole field is also reducible to 
the two-component radial equation l|5(J|) . with k — (see [28| for technical details). 

We now show that Eq. I|5U|) is factorizable when V = in the sector j = \q\ — 1/2 (k = 0). In fact, if the wave 
function is transformed to 

n ^-rt ({) , T = \,( \ ]) , ,54, 



9 J V^+y V 9 J ' V2\- 11 

then the Dirac equation (|50|l changes to the following factorized form 

d 



— ~r C/(r)J V- - ^+ , 

~ + I/(r) ) ^+ = EiP_ . (55) 
ar / 

Here U{r) plays the role of the superpotential W(r), and the energy parameter is e = E 2 . The form of U{r) admitting 
exact/quasi-exact solutions can then be constructed according to the prescribed procedure. 

A. Exactly solvable cases 

From the Table 4.1 of it is found that there are four types of exactly solvable potentials U(r) for this system, 
namely, (i) 3D oscillator-like ; (ii) Coulomb-like; (iii) Eckart potential; and (iv) generalized Poschl- Teller. 

B. Quasi-exactly solvable cases 

Three types of QES potentials U{r) can be constructed, namely, Class VII, VIII, and IX in Q. We outline the 
construction of Class VII below. 

The general potential in Class VII has the form 

V N (r) = a 2 r 6 + 2abr 4 + [b 2 - a (4JV + 2 7 + 3)] r 2 + 7 (7 - 1) r~ 2 - b (2j + 1) , (56) 

where a, b and 7 are constants. The gauge function is 

g{r) = ^r 4 + ^r 2 - 7lnr . (57) 

We must have a, 7 > to ensure normalizability of the wave function. The potential 

U (r) = g'{r) = ar 3 + br- - (58) 

r 

admits a QES ground state with energy E = 0, and the corresponding wave function with components tp- oc 
exp(-5oM) and ip+ = 0- 

To determine potentials admitting QES potentials Vn with higher weights, we need to obtain the Bethe ansatz 
equations for <fi. Letting z = h(r) = r 2 , Eq. becomes 



-4^—5- + Uaz 2 + 46z - 2 (2 7 + 1)) — - UaNz + e) 
dz z dz 



(z) = 0. (59) 



For N = 0, the value of the e is e = 0. For higher N > and </>(r) = riili( z — z i)i the potential U^{r) is obtained 
from Eqs. (|Tg|l : 

Mr) = W-E^. (60) 



For the present case, the roots Zi's are found from the Bethe ansatz equations 

_ Z; 



2az 2 + 2bz l -{2 1 + l)- s T ^— = , i = l,...,N, (61) 

^— ' Zi - Zl 



and e in terms of the roots Zj's is 



N 1 

e = 2(27 + 1)2- ■ ( 62 ) 



(=i 



. (63) 



2a 



e ± = 2(6± V^ 2 + 2a(27+l)) . (64) 



For iV = 1 the roots z\ are 



and the values of e are 



For a > 0, the root = — |z 1 | < gives the ground state. With this root, one gets the potential 

Ux{r) =ar 3 + br-— — t — --2. (65) 
r 2 + \z{ | r 

QES potentials for higher degree N can be constructed in the same manner. 

VI. SUMMARY 

The exact /quasi-exact solvability of the Dirac equation with a Lorentz scalar potential is examined in this paper. 
Possible forms of the potentials allowing factorization of the Dirac equation are identified in two, three, and four 
dimensions. Based on such factorizability, we have classified all exactly solvable potentials, and QES potentials based 
on sl(2) Lie-algebra according to Turbiner's calssification. 

In this work and in [ITIIT^ |. we have classified exactly solvable and QES potentials in the Pauli equation and the Dirac 
equation coupled minimally to a stationary vector potential, the Dirac-Pauli equation (equivalent to generalized Dirac 
oscillator) coupled non-minimally to an external electric fields, and Dirac equation with a Lorentz scalar potential. 
Nonetheless, we should mention that we have by no means exhausted all possibilities of exact/quasi-exact potentials. 
Our classification and construction are based on the factorizability, or supersymmetric structure, of the system, with 
sl(2) as the underlying symmetry. There are QES systems which are not factorizable, such as those considered in 
[ii. Il3l ITrj Il7 | . and systems which are not related to the Lie-algebra sl(2), such as that discussed in It is an 

interesting and challenging task to develop new methods to classify and construct QES potentials in multi-component 
wave equations with any Lie-algebraic structures. 
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